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We study the conditions under which thermal fluctuations generated in the contracting phase of a 
non-singular bouncing cosmology can lead to a scale-invariant spectrum of cosmological fluctuations 
at late times in the expanding phase. We consider point particle gases, holographic gases and string 
gases. In the models thus identified, we also study the thermal non-Gaussianities of the resulting 
distribution of inhomogeneities. For regular point particle radiation, we find that the background 
must have an equation of state w = 7/3 in order to obtain a scale-invariant spectrum, and that 
the non-Gaussianities are suppressed on scales larger than the thermal wavelength. For Gibbons- 
Hawking radiation, we find that a matter-dominated background yields scale-invariance, and that 
the non-Gaussianities are large. String gases are also briefly considered. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

In recent years, there has been a lot of interest in non- 
singular bouncing cosmologies (see e.g [l| for a recent 
review with an extensive list of references). Such cos- 
mologies may be desirable since they resolve the singu- 
larity problem of the inflationary scenario, the current 
paradigm of early universe cosmology. If inflation is re- 
alized by making use of the potential energy of a scalar 
matter field while treating space-tiime dynamics using 
the Einstein action, then an initial cosmological singu- 
larity is unavoidable Q. 

While it is possible that that in the context of a non- 
singular bouncing cosmology a period of inflationary ex- 
pansion [3] is realized after the bounce (see e.g. [J] for 
an explicit model) and the cosmological perturbations 
observed today are generated as quantum vacuum fluc- 
tuations in the inflationary phase [5], it is also possible 
to obtain fluctuations without requiring an inflationary 
phase after the bounce. The reason is that all scales ob- 
served today were at one point inside the Hubble radius 
at sufficiently early times during the contracting phase. 
This is illustrated in Figure 1, which shows the space- 
time plot in a non-singular bouncing cosmology [431 ] . 

In the context of studies of bouncing cosmologies it 
has been realized that fluctuations which are generated 
as quantum vacuum perturbations and exit the Hubble 
radius during a matter-dominated contracting phase lead 
to a scale- invariant spectrum of cosmological fluctuations 
today [ll|, [13, Ell ( see a l so 0X an d thus yield an alter- 
native to cosmological inflation for explaining the current 
observational data. 



However, in a bouncing cosmology it is not manifest 
that perturbations arise as quantum vacuum fluctua- 
tions. In inflationary cosmology it can be argued that 
the exponential expansion of space during the inflation- 
ary phase red-shifts all classical matter initially present 
and leaves behind a matter vacuum. However, if the uni- 
verse starts out large and cold in a contracting phase, 
there does not seem to be a reason to single out vacuum 
over thermal initial conditions for the fluctuations. 

Thermal fluctuations as the origin of structure in the 
universe were considered in the context of a forever ex- 
panding cosmology, but it was concluded that it was not 
possible to obtain a scale-invariant spectrum of cosmo- 
logical perturbations [HI]. However, as realized in 0], 
thermal fluctuations of strings can give rise to a scale- 
invariant spectrum if we abandon the assumption that 
the universe is forever expanding. Specifically, in [7J it 
was shown that thermal string fluctuations in a quasi- 
static early string phase yield a scale- invariant spectrum. 
Hence, there are good reasons to expect that in a bounc- 
ing cosmology it might be possible for thermal fluctu- 
ations to generate a spectrum whose shape is in good 
agreement with observations. 

In this paper we will study the conditions under which 
thermal initial conditions for perturbations can lead to 
a scale-invariant spectrum of cosmological fluctuations 



after the bounce [44j. We find various possibilities, de- 



* caiyf@ihep.ac.cn 
txuewei@physics.mcgill.ca 
trhb@physics.mcgill.ca 
§ xmzhang@ihep.ac.cn 



pending on what kind of thermal fluctuations we consider 
If the fluctuations are generated by normal particle 
radiation (with the usual equation of state w r = 1/3, 
where w is the ratio w = p/p of pressure p to energy 
density p), then the required equation of state of the 
background is w — 7/3. For radiation in holographic 
cosmology (ThI . we find that a background equation of 
state of w = will lead to a scale-invariant spectrum. 
String gases require a quasi-static early phase. 

As realized [l{| in the study of perturbations in the 
"matter bounce" , in which vacuum fluctuations exit the 
Hubble radius during a matter-dominated phase of con- 
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traction, inhomogeneities generated during a contracting 
phase leave a distinctive imprint on the magnitude and 
shape of the non-Gaussianities manifest in the three point 
function. We will study these non-Gaussianities for ther- 
mal initial conditions. 

The outline of this paper is as follows: In the fol- 
lowing section we study the evolution of fluctuations in 
a bouncing universe. In Section 3, we consider ther- 
mal initial conditions for the perturbations and ask un- 
der which conditions a scale-invariant spectrum after the 
bounce results. In Section 4 we then estimate the non- 
Gaussianities in the resulting models. 



II. FLUCTUATIONS IN A BOUNCING 
UNIVERSE 

It is useful to take a first glance at the route of cos- 
mological perturbations in a bouncing cosmology. As is 
depicted in Figure [TJ the wavelength of fluctuations be- 
comes larger than the Hubble radius in the contracting 
phase, and reenters in the expanding phase. All that is 
required for this space-time sketch to apply is that the 
equation of state parameter w of the background universe 
be larger than —1/3 at all times except possibly around 
the bounce point. This requirement is automatically sat- 
isfied for most forms of matter (the exception being scalar 
field models which lead to inflation) . In Figure [1] we plot 
the evolution of the physical length corresponding to a 
fixed comoving scale. This scale is the wavelength of the 
fluctuation mode k which we will follow later. 
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case, the linearized fluctuations about a Friedmann- 
Robertson- Walker background metric in longitudinal 
gauge can be expressed as (see e.g. [2(| for a review 
of the theory of cosmological perturbations) 



a(v) 2 [{l + 2<£)dr] 2 - (1 - 2<$>)dx 2 



(1) 



where $(a;, ?y) is the generalized Newtonian gravitational 
potential which characterizes the metric fluctuations, r\ 
is conformal time and a{rf) is the background scale factor 

In Fourier space, the fluctuation variable $ satisfies 
the following second order differential equation 



2aH& k + (c 2 s k 2 - 2eH 2 + 2oH 2 )<$> k = , (2) 



where TL = a' /a is the comoving Hubble parameter, and 
the prime denotes the derivative with respective to 77. 
The parameter c s is the sound speed, which we take to be 
a free parameter. Moreover, we have defined two useful 
parameters 
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and 



(3) 
(4) 



which characterize the background evolution. For a con- 
stant equation of state w, we have 



(5) 



and thus the perturbation equation can be greatly sim- 
plified. 

We will consider an equation of state for which 

the background scale factor evolves as a power of time. 
If w > —1/3, an equation of state which does not give 
accelerated expansion, we can choose our origin of the 
time axis such that t = corresponds to the bounce 
point. More specifically, we shall take t = to be the 
value of t for which the Big Crunch singularity would 
occur in the absence of the terms in the action which lead 
to a non-singular bounce. In this case, we can also choose 
the origin of the conformal time coordinate such that t — 
corresponds to 77 = 0. With this choice of coordinates, 
then in the contracting phase, the perturbation equation 
@ takes the form 



l + 2v 



c 2 s k 2 $ k = , 



(6) 



FIG. 1: A sketch of the evolution of scales in a bouncing uni- 
verse. The horizontal axis is the physical spatial coordinate, 
the vertical axis is time. Plotted are the comoving Hubble ra- 
dius |7^| _1 and the wavelength A of a fluctuation mode with 
fixed comoving wavenumber k. 

In this paper we focus on adiabatic fluctuations and 
consider matter without anisotropic stress. In this 



where 



5 + 3w 
2(1 + 3w) 



(7) 



Since © is a second order differential equation, there 
are two linearly independent solutions. On super-Hubble 
scales, one mode is constant (the "D-mode"), the other is 
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growing in a contracting phase (and decreasing in an ex- 
panding phase) . We call the second mode the "S-mode" . 
The general solution on super-Hubble scales is a linear 
combination of the two modes and is hence given by 



_S- 

(V) 



2// 



(8) 



where Z?_ and SL are the mode coefficients. If the equa- 
tion of state is the same in the expanding phase after the 
non-singular bounce, then the solution in the expanding 
period can be written as 



S+ 

{rj) 2u 



(9) 



with new mode coefficients D + and S + . In a non- 
singular bounce, the fluctuations can be smoothly 
evolved through the bounce. Thus, the mixing ma- 
trix which relates the mode coefficients in the expanding 
phase with those in the contracting phase can be calcu- 
lated m 

The key question which arises when studying the 
transfer of fluctuations through a non-singular bounce is 
whether the spectrum of the dominant mode in the con- 
tracting phase, the S- mode, couples to the D + mode, 
the dominant mode in the expanding phase. This is- 
sue was initially studied by replacing the bounce region 
with a matching surface across which the perturbations 
are connected making use of the Hwang- Vishniac [23| 
(Deruelle-Mukhanov [23]) matching conditions. It was 
found [IE [IE [13] that the coupling is suppressed by a 
factor of k 2 on large wavelengths, i.e. 



Da 



k 



0(1)D_+0(1)(-YS- 



(10) 



where k* is a normalization scale which is set by the mi- 
crophysics, i.e. is in the ultraviolet range. However, as 
pointed out in [28], it is not valid to apply the matching 
conditions to fluctuations at an interface between a con- 
tracting and an expanding phase because the background 
does not satisfy the matching conditions. If the bounce 
is non-singular, however, one does not need to make 
use of matching conditions: the perturbations can be 
evolved continuously from the contracting to the expand- 
ing phase. This was done in the context of a non-singular 
regularization of Pre-Big-Bang cosmology in [2^] , and in 
the context of a non-singular Ekpyrotic model in [30] ]. 
More recently, such calculations were carried out for a 
non-singular mirage cosmolo gy 13 II] . in the non-singular 
higher derivative bounce of [l6| in [32j], for a quintom 
bounce model i n [331) . and more specifically in the Lee- 
Wick bounce in [22[ . 

For a non-singular bounce, one can follow the fluc- 
tuations both numerically and analytically. To obtain 
a good analytical approximation, one divides the back- 
ground time into three intervals - the contracting phase, 
the bounce phase where the Hubble expansion rate can be 
modelled as H (t) — at (where a is a constant whose value 



is set by the new physics which determines the bounce), 
and the post-bounce expanding phase. The duration of 
the bounce phase is set by the scale determining the new 
physics which regulates the bounce. The result of the 
works quoted above is that on length scales larger than 
the duration of the bounce, the mode mixing occurs as 
given by (|10p . Thus, the contribution of the dominant 
mode in the contracting phase is suppressed in the dom- 
inant expanding phase mode function by a factor of k 2 . 
On the other hand, on scales short compared the bounce 
time, there is no suppression of the coupling. 

In the following, we will, without much loss of gen- 
erality, assume that the bounce is short from the point 
of view of cosmological scales of interest. Thus, we will 
use the coupling given by (flO]) . In [22j], we studied the 
evolution of perturbations which start out as quantum 
vacuum fluctuations and discovered that it is precisely 
quantum vacuum fluctuations which exit the Hubble ra- 
dius is a matter-dominated contracting phase which are 
scale-invariant after the bounce. However, in the frame- 
work of a cosmological model which starts out large and 
cold, there is no particular reason to focus on quantum 
vacuum initial fluctuations. It is rather reasonable to 
consider initial thermal perturbations. In this paper we 
will study under which conditions on the background cos- 
mology one obtains a scale-invariant spectrum at later 
times starting from thermal initial conditions. 

Before starting the analysis, we remind the reader of 
a useful formula for the time of Hubble radius crossing. 
Given an equation of state parameter w, the scale factor 
evolves as 



i(t) ~ t p 



with p = 



3(1 



(11) 



and yields the following relation for the conformal time 
V 



n 



t i- P 



The condition for Hubble radius crossing is 
k ~ aH 



(12) 



(13) 



for a perturbation mode with fixed comoving wavenum- 
ber k. This yields 



and thus 



rj H (k) ~ k 1 , 



t H {k) ~ k 



(14) 



(15) 



where the subscript "U" denotes the moment of Hubble 
radius crossing. 



III. THERMAL FLUCTUATIONS 

The method of calculating the spectrum of cosmologi- 
cal perturbations at late times is the following. First, we 
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compute the matter fluctuations on sub-Hubble scales in 
the contracting phase. Next, for any scale k, we compute 
the induced metric fluctuations at the time tn{k) when 
the scale exits the Hubble radius during the contract- 
ing phase. In the third step, the metric fluctuations are 
evolved on super-Hubble scales making use of the evolu- 
tion equations for perturbations discussed in the previous 
section. This is the standard way of following the gen- 
eration and evolution of cosmological perturbations, as 
app lied to inflationary cosmology in early works (see e.g. 
[Ij] and [35[) and to string gas cosmology @, S 

hj. The 

method of calculation reflects the fact that metric fluctu- 
ations are sub-dominant on sub-Hubble scales, but that 
on super-Hubble scales the matter fluctuations freeze out 
and the evolution of the perturbations is driven by the 
metric. 

The key constraint equation which relates matter and 
metric fluctuations is the time-time component of the 
perturbed Einstein equation 

-m{H® + &) +V 2 $ = 4nGa 2 5p , (16) 

where Sp is the fluctuation of the energy density. Note 
that all the three terms on the left-hand-side of the above 
equation are of the same order of magnitude at the Hub- 
ble radius crossing time. Therefore, up to a constant of 
order O(l), the power spectrum of the metric perturba- 
tions is given by, 



P*(Je) 



1 



12tt 



;k 3 \mf 



(17) 



1 _ fc3 < Sp(k) 2 > 



4M p 4 (H(t H (k)) 4 



where (in our case) the pointed brackets denote en- 
semble averaging in thermal equilibrium. Making use 
of the Hubble radius crossing condition H(tff(k)) — 
a ■ {tn(k))k and replacing the power spectrum of the 
Fourier space energy density correlation function by the 
position space correlation function we obtain: 

P*(k)(t H (k)) = ^lj*- 4 o 4 (^(fc)) < Sp 2 > \ m , 

(18) 

where < Sp 2 > \ntk) (t° be evaluated at Hubble ra- 
dius crossing) is the position space energy density fluc- 
tuation correlation function is a sphere of radius R(k), 
where R(k) is the physical length corresponding to the 
co-moving momentum scale k. 

In a system which is in thermal equilibrium, the cor- 
relation function of the energy density is given by 



< Sp 2 > | fi(fe) = C V (R) 



11 
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(19) 



where Cv(R) is the heat capacity in a sphere of radius 
R and is defined in terms of the expectation value of the 
internal energy as 



C V (R) = §f<E> 



(20) 



A. Thermal Particle Fluctuations 

In this subsection we consider fluctuations in a gas of 
point particles with an arbitrary equation-of-state w r . In 
this case, from the stress-energy conservation equation 
it follows that the energy density and the temperature 
change as a function of the scale factor as follows: 



p r ~ a 
T ~ a~ 3w 



-3(l+3uv) 



which yields 

Pr ± j 

and so we get the heat capacity 



Cv{R) = R 3 % ~ R 3 T^ 



dT 



(21) 
(22) 



(23) 



(24) 



Inserting (JT^ and (J24|) into (TTT]) and applying R ~ 
1/H, we obtain the power spectrum for metric perturba- 
tions at the Hubble radius crossing time tjj(fc): 



l-3p(l + 2i+y) 



(25) 



P 9 (k) ~ T H ^H~ H \ k) ~ k~ 
where we have used the relations 

T H (k) - a- 3w - - k^f , and (26) 



H t H (k) ^ k^-p 



(27) 



Since the power spectrum of the constant mode Pd (k) 
is the same as that of $ at Hubble radius crossing, it 
scales as 



P D {k) 



no = 



k riD 

1 - 3p(l + 2w r ) 



(28) 



The power spectrum of the growing mode Ps(k) is the 
spectrum of $ at Hubble radius crossing modulated by 
the factor r]H(k) 4 " . This yields 



Ps(k) 



n s = 



P D (k)k- 4v 
3-p(l 



k ns 



(29) 



6w r 



p-1 



There are two possibilities to obtain a scale-invariant 
spectrum after the bounce. The first is if the D-mode in 
the contracting phase is scale-invariant (recall that the 
contribution of the contracting phase D-mode to the ex- 
panding phase D-mode is not suppressed). The second is 
that the S-modc in the contracting phase has a spectrum 
proportional to fc~ 4 . Since the contribution of the con- 
tracting phase S-mode to the expanding phase D-mode 
is suppressed by k A in the power spectrum (see (1101) , this 
then yields scale-invariance of the dominant mode in the 
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expanding phase. Thus, the two possibilities to obtain a 
scale-invariant spectrum are 



n D = , n s + 4 > ; 
or no > , ns + 4 = 



(30) 
(31) 



In order for the D-mode to provide a scale-invariant 
spectrum, we require 



w r = -(w - 1) 



(32) 



This is a solution provided that the S-mode yields a blue 
spectrum which requires 



1 



1/3 



> 



(33) 



which is satisfied if w > 1. Note that for such a back- 
ground, the anisotropic stress grows less fast than the 
background energy density. Thus, this background is sta- 
ble towards anisotropic stress fluctuations. Note also that 
if we were to demand that the equation of state of the 
background and of the fluctuations is the same, we would 
obtain w = —1/3 which is not a physical solution since 
it is the borderline solution between an acelerating and a 
decelerating background, and no scales exit the Hubble 
radius during the contracting phase. 

For the growing mode (S-mode) in the contracting 
phase to provide a scale-invariant spectrum in the ex- 
panding phase, the condition on the two equations of 
state is 



w r = -(w-1) , 



and the dominance of the S-mode requires 



(34) 



(35) 



In the special case that the radiation is normal radia- 
tion with w r = |, the only possibility to obtain a scale- 
invariant spectrum is via the D-mode in the contracting 
phase. This requires 



IV 



(36) 



B. Fluctuations in Holographic Cosmology 

Gibbons-Hawking radiation (36|, originally discovered 
in studies of thermodynamics in de-Sitter space, has re- 
cently been studied extensively in the context of devel- 
opments in string theory, in particular in light of the role 
of holography in string theory [l8[ . The key point is that 
quantities such as the entropy and the energy do not scale 
extensively with the size of the volume, but increase as 
the area [18|]. Thus, the average energy is 



< E > = TR 2 M 2 



(37) 



which gives a special heat capacity 
Cv(R) ~ R 2 M 2 



(38) 



If we combine Eqs. flTjl. (fT9|) and l(38]). and use the def- 
inition of the Gibbons-Hawking temperature associated 
with the instantaneous Hubble radius 



T 



1 

R 



(39) 



where R is taken to be the Hubble radius R ~ l/-ff, then 
the power spectrum for $ at the Hubble radius crossing 
is given by 



P*(Mff(fc)) ~ ( 



T H (k) 



(40) 



H(t H (k)) 2 
M„ 1 



fc2/(X-p) 



Thus, a scale-invariant spectrum from the D-mode can 
be achieved in the limit p — > oo, which corresponds to 
an inflationary contraction. However, in this case scales 
are not exiting the Hubble radius during the contracting 
phase, and thus it does not make sense to consider ther- 
mal fluctuations in this context since thermal equilibrium 
cannot be established on super-Hubble scales. 
The contribution of the S-mode scales as 



;+2 



(41) 



Inserting the relation for v in terms of p we find that 
the condition for scale-invariance is 



= 0. 



(42) 



As a consistency check, we note that, as follows from 
PO")) . for p = 2/3 the D-mode has a red spectrum. Hence, 
this is indeed a background equation of state for which 
thermal fluctuations of a holographic gas yields a scale- 
invariant spectrum. 



C. Thermal Fluctuations of a String Gas 

For completeness, we will also add an analysis for ther- 
mal string gas fluctuations. If space is compact, then the 
specific heat capacity of a gas of closed strings also sat- 
isfies the holographic scaling 



C V (R) 



R 2 M 2 pl 



(43) 



Inserting this relation into (fT8|) . ([T9| and (|20|) we find 

P*(h,t„(k)) ~ a 4 (t H (k))T 2 (t H (k)). (44) 

In the case of string gas cosmology jf| 0, 0, E3] , fluctu- 
ations exit the Hubble radius at the end of a quasi-static 
Hagedorn phase. Thus, both a(tjj(k)) and T(tji(k)) are 
almost independent of k, and a scale-invariant spectrum 
(with a slight red-tilt) results. 
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If we forget about the Hagedorn background (which 
cannot be described in terms of the Einstein or dila- 
ton gravity background equations), and simply couple a 
string gas to the non-singular background geometry dis- 
cussed in this paper, and use T(t) ~ a(i) _1 , then the 
D-mode of $ yields a scale- invariant spectrum if p = 0. 
This is consistent with the Hagedorn phase of string gas 
cosmology. We also find that the S-mode yields a scale- 
invariant contribution to the post-bounce spectrum of $ 
if p = 1/4, but for this value of p, the contribution from 
the D-mode has a red tilt and hence dominates. 



IV. NON-GAUSSIANITIES FROM THERMAL 
FLUCTUATIONS 

Due to the non-linearities in the theory, the fluctua- 
tions are not perfectly Gaussian. A lot of recent interest 
has focused on calculating the non-Gaussianities as man- 
ifested in the three-point function (see e.g. [38] for a 
review). In single field slow- roll inflation models the am- 
plitude of the predicted non-Gaussianities is suppressed 
by the slow-roll parameter. In models with a contracting 
phase, however, the induced non-Gaussianities are typi- 
cally much larger (see e.g. [39| for studies in the context 
of the Ekpyrotic scenario). 

As studied recently in [l9[, the non-Gaussianities as 
measured by the three-point function are of order one 
in the non-singular matter bounce scenario in which the 
fluctuations are of quantum vacuum origin. The key 
facts that lead to this result are firstly that there is 
no suppression of the non-Gaussianities by slow-roll pa- 
rameters, and secondly that fluctuations grow on super- 
Hubble scales in the contracting phase, which results in 
a larger time interval determining the amplitude of the 
effects and in different terms dominating the shape of the 
three point function. 

In this section we calculate the non-Gaussianity esti- 
mator /jvl in bouncing cosmologies with thermal fluc- 
tuations. We will use the formalism to compute non- 
Gaussianities of thermal origin which has been devel- 
oped in [40l | and which was applied to estimate the non- 
Gaussianities in strin g ga s cosmology in [4l[ and to holo- 
graphic cosmology in [42| • We will calculate the following 
non-Gaussianity estimator for fluctuations on a scale k: 



/jvi(fc) 



<ct> 



(45) 



18Vfc3 < Q > A 

where £ is the curvature fluctuation in co-moving gauge. 
Thus if we obtain the two-point and three-point corre- 
lators of metric perturbations originating from thermal 
fluctuations, the non-Gaussianity parameter can be cal- 
culated |47j . 

The expression for the two-point correlation function 
was given in (|19[) . The three-point correlation function 
of thermal fluctuations in an equilibrium ensemble is 



< Sp 3 > = - 



1 dHnZ = T^d_ (r 2] 
R 9 dp 3 R 9 dT R ' 



(46) 



where Z is the partition function and /3 is the inverse 
temperature. 

A. Normal Radiation 

For normal radiation with equation of state w r = 1/3, 
we need the background equation of state of the universe 
in the contracting phase to be w — 7/3 in order to obtain 
a scale-invariant spectrum of fluctuations, as studied in 
Section 3. In this case, the heat capacity can be expressed 
as 



C V (R) = c v R 3 T 3 



(47) 



where c v is determined by the background initial con- 
ditions and here is treated as a constant of order O(l). 
From this, we obtain the following expressions for the 
two-point and three-point correlation functions of the 
density perturbations 



<6p 2 > = 
<5p 3 > = 



„T 5 



R 3 1 
5c v T 6 
~R^~ 



From Eq. (|17[) we find the relation 



Sp 



2M£H 2 



(48) 
(49) 

(50) 



which is valid for any thermal system. In Section 3 
we have seen that for the background considered here, 
the dominant mode after the bounce is seeded by the 
constant mode (the D-mode) in the contracting phase. 
Hence, the relation between $ and £ is given by 



c 



3u> 



3 + 3w 



4>; 



G 



4>; 



(51) 



Combining the above equations, and going from 
Fourier space to position space via 

Sp = 2-^ 2 n- 1 5p(k) , (52) 

where the left-hand side represents the root mean square 
density fluctuation corresponding to the co-moving scale 
k, finally the non-Gaussianianity estimator takes the 
form 



25 



Tnl 



< Sp 3 > 



(53) 



where the right hand side is evaluated at the Hubble 
radius crossing time. Inserting the above expressions for 
the density two and three point functions we get 



125 



TNL 



54V2^c,T 4 (^(fc)) 



(54) 
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Noticing that for w r — 1/3 we have T oc a 1 , inserting 
H{tn{k)) ~ (&)) and making use of (fT5|) we obtain 



iVL 



fc 



(55) 



and the coefficient is such that for scales exiting the Hub- 
ble radius right before the bounce (when all quantities on 
the right hand side of (I54[) are of the order of the Planck 
mass - assuming that the energy scale at the bounce is 
given by the Planck mass) the amplitude of /jvz, is of or- 
der 1. Inserting the value p = 1/5 obtained for w = 7/3 
we finally obtain 



f NL (k) 



0(l)(f- f/\ 

Kb 



(56) 



where ks is the value of k for which the wavelength exits 
the Hubble radius immediately before the bounce. 

The above analysis shows that the non-Gaussianities 
in a bouncing cosmology in which the fluctuations are 
seeded by particle thermodynamic perturbations are 
Poisson suppressed on wavelengths larger than those ex- 
iting the Hubble radius immediately before the bounce. 
Thus, the non-Gaussianities are very different from what 
is obtained [191 ] in the case of a matter bounce with quan- 
tum vacuum initial perturbations. Besides the Poisson 
suppression of thermal non-Gaussianities (by the central 
limit theorem, on large scales the thermal fluctuations 
must approach a Gaussian), an important reason for the 
difference is that in the case of the "matter bounce" of 
the S-mode is responsible for the final fluctuations. 
The S-mode grows on scales larger than the Hubble ra- 
dius in the contracting phase, thus leading to an enhance- 
ment of the non-Gaussianities. 



B. An Estimate of /nl for Gibbons-Hawking 
Radiation 



In the case of a holographic gas we expect the resulting 
non-Gaussianities to be much larger since for each scale 
the thermal correlation length is taken to be equal to the 
Hubble radius at the time that the scale exits the Hubble 
radius. Hence, we do not expect a Poisson suppression 
factor, and we expect non-Gaussianities to be of order 1. 
This expectation is verified by an explicit computation. 

The starting point is the expression Cy{R) — c v M 2 R 2 
for the heat capacity, which leads to the following expres- 
sion for the two-point and three-point correlation func- 
tions of energy density perturbations 



<6p 2 > 



< Sp 3 > 



c v M 2 T 2 

i? 4 ' 
2c v M 2 T 3 

R 7 



Inserting these expressions into (|53p we obtain 



TNL 



25 1 RH 2 
54V2tt c v T 



(57) 
(58) 

(59) 



where all quantities on the right hand side are to be eval- 
uated at Hubble radius crossing. 

Making use of the Gibbons-Hawking relation T ~ 
1/R ~ H, we then have 



f ~ 1 

JNL — — 

C v 



(60) 



which shows the non-Gaussianity is of order 0(1) and 
scale independent. 

Note that the sign of Jnl is positive in the case of 
thermal fluctuations considered here. This is different 
from the negative sign obtained in the matter bounce 
scenario with vacuum initial conditions [l||. 



V. CONCLUSIONS 

In this paper we have studied the possibility of obtain- 
ing a scale-invariant spectrum of fluctuations from ther- 
mal initial conditions in the context of a non-singular 
bouncing cosmology. We classified the conditions on the 
equation of state w r of the thermal radiation and the 
equation of state w of the background which yield a scale- 
invariant spectrum. 

In the case of regular particle radiation with w r = 1/3 
we find the condition w — 7/3 for the background. Such 
as equation of state can be realized by supposing that the 
background is determined by a scalar field with a nega- 
tive exponential potential, similar to what is assumed in 
Ekpyrotic cosmology [37| • In this case, the fluctuations in 
radiation would be entropy fluctuations which would then 
seed adiabatic fluctuations with the same spectral index 
(note that the primordial adiabatic fluctuations induced 
by the scalar field are blue and hence will not dominate 
for long wavelength modes [25ll26T|). 

In the case of a Gibbons-Hawking radiation we find 
that an equation of state w = (a matter bounce) yields 
a scale-invariant spectrum of cosmological perturbations. 

Finally, to obtain a scale-invariant spectrum of cosmo- 
logical perturbations with a string gas requires a quasi- 
static early phase. This matches with what is usually 
assumed in string gas cosmology fToj ] . 

We have also considered the non-Gaussianities of ther- 
mal fluctuations. In the case of thermal particle fluc- 
tuations we find that the non-Gaussianities are Poisson- 
suppressed on large scales. They are of the order 1 only 
on microscopic length scales, scales for which the thermal 
correlation length at the time of Hubble radius crossing 
is comparable with the Hubble radius itself. However, 
in the case of radiation satisfying the Gibbons-Hawking 
distribution, the non-Gaussianity is of order 1 /c v (where 
c v is the constant which determines the specific heat ca- 
pacity) and approximately scale independent. 
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